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The problems associated with practical implementation of the scheme proposed for preparation 
of arbitrary states of polarization ququarts based on biphotons are discussed. The influence of 
frequency dispersion effects are considered, and the necessity of group velocities dispersion com- 
pensation in the frequency non-degenerate case even for continuous pumping is demonstrated. A 
method for this compensation is proposed and implemented experimentally. Physical restrictions 
on the quality of prepared two-photon states are revealed. 
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So-called high-dimensional quantum systems or systems with a dimension of Hilbert space D > 3 have 
attracted increasing interest. Primarily, this interest stems from the possibility of using such systems in 
^C) \ various quantum information protocols. The application of these systems, which are frequently referred to 

as audits (quantum dits) , as information carriers offers a number of advantages over encoding with qubits 
[l|-|4|]. The second important field of application of qudits is verification of basic principles of quantum 
theory, in particular, violation of Bell-type inequalities [f|. Incidentally, it turns out that the use of high- 
dimensional systems in a number of cases makes it possible to increase the quantitative gap that arises 
in the description of composite systems using the language of quantum or classical correlations. In turn, 
' this allows one to reduce the requirements imposed on quantum efficiency of photo-detectors and other 

experimental equipment Thus, the development of methods used for preparing and characterizing 
different states of qudits is an important problem giving rise to a new field of quantum information, i.e., 
so-called quantum state engineering. 

Polarization degrees of freedom of photons are especially convenient tools for experimental preparation 
^ ' of qudit states, primarily owing to simplicity of performing polarization transformations. Specifically, 

C*~) , in order to control the polarization state, it is sufficient to use linear optical elements, like phase plates 

and beam-splitters. Moreover, the developed statistical methods of polarization tomography have made 
it possible to completely reconstruct the state vector of the initial quantum system as a result of the 
relatively simple procedure 0, § . 

A polarization ququart is defined to be a pure two-photon state that can be written in the form 

O: |^>=Ci|fl- 1 )|fr 2 )+C3|Fi)|y a >+C3|V r i>|^2>+C4|Vi)|^2>, (i) 
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where the ket vectors \H) and \V) stand for the basis single-photon polarization states of the field mode, 
and indices 1 and 2 can refer to either frequency (longitudinal) or spatial (transverse) modes. It is com- 
mon practice to distinguish spatial polarization ququarts, for which the field state is quasi-degenerate in 
frequency but there are two spatial modes (see, for example, Q), and frequency polarization ququarts, 
when the spatial modes are degenerate but the frequency spectrum of the field contains two components 
Q . The first experiments on the preparation of ent ang led pairs of polarization qubits or entangled states 
of ququarts were performed in the mid-90s. In fiol \1M pairs of photons generated in collinear frequency- 
degenerate spontaneous parametric downconversion (SPDC) with type-II phase-matching were trans- 
formed into polarization-momentum entangled states with a beam splitter and subsequent postselection. 
Shortly thereafter, Kwiat et al. [H] proposed a method that does not require postselection and is based 
on the type-II non-collinear phase matching. There are well-known schemes employed for preparation 
of entangled states with the use of two crystals with type-I phase-matching in non-collinear frequency 
degenerate [Til and collinear but frequency non-degenerate [HI regimes. The case of pulsed pump is 
discussed in |l6l [l7| . One of the advantages of these schemes is the possibility of preparing biphotons in 
non-maximally entangled states [l8j . Of special interest, however, are the methods providing the prepa- 
ration of polarization ququarts in an arbitrary state by means of unitary (lossless) transformations, when 
all components of the polarization state of a biphoton pair can be completely controlled. The design 
of this universal source of ququarts is of both fundamental and practical interest for solving problems 
of quantum information and quantum communication. The experimental scheme that ensures complete 
control of a polarization ququart was proposed in a recent paper [T^ | . Being essentially an interferomet- 
ric scheme it provides the stability of parameters, which is not inherent in conventional interferometric 
schemes but is well known in classical polarization optics. To the best of our knowledge, this scheme 
has not been fully implemented so far. However, specific families of states of frequency non-degenerate 



ququarts, of course, have already been prepared experimentally (see, for example 



quency non-< 

, SEE Hi). 
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In this work, we study physical restrictions imposed on the quality of states prepared according to 
the aforementioned scheme and discuss the methods proposed for their elimination. It has turned out 
that the use of the frequency non-degenerate regime and, correspondingly, the polarization and frequency 
entanglement in the scheme with two consequent crystals involves specific experimental difficulties due 
to the frequency dispersion in crystals. However, this scheme is very convenient from the viewpoint of 
applications in quantum information protocols and has been frequently employed in their implementation. 
In this paper, we demonstrate that, in order to achieve high quality of prepared states, it is necessary 
to use compensators for group velocity dispersion, which were earlier treated as an attribute of schemes 
with type-II phase matching [lOL HP . l2l| or pulse pumping [l6l . [l7| only. 

II. SCHEME FOR PREPARATION OF ARBITRARY POLARIZATION STATES OF 

QUQUARTS 

Let us consider the scheme proposed in [l9j for the preparation of arbitrary states of polarization 
ququarts. We start with the formal mathematical aspect of the problem. It is known that an arbitrary 
pure entangled state of a qubit pair or a ququart can be represented in the form (cf. expression (J]}) 

|V>) = y/Jl \Ai) \A 2 ) + ^j—JI |B0 \B 2 ) , (2) 

where \Aj) and \Bj) are the basis vectors in single-qubit state spaces. This expression is known as the 
Schmidt decomposition. The coefficients /x, (1 — fj.) are eigenvalues of the single-particle density matrices 
of each qubit (which, as is known, coincide with each other), and the vectors \Aj) and \Bj) form an 
orthogonal basis in which they are diagonal. Therefore, in order to prepare an arbitrary ququart state, 
one should be able to experimentally control the coefficients in the Schmidt decomposition and to perform 
switching between the bases. Note that the former operation changes the degree of entanglement of the 
state and affects the degrees of freedom of both qubits. In this sense, such an operation can be termed 
"nonlocal" [.30] . However, the transformation of the basis vectors requires only "local" operations, i.e., 
operations that act on each qubit separately. It should also be noted that the degree of entanglement of 
a_^jubit pair can be defined in any reasonable manner; for example, it is convenient to use concurrence 

The idea of possible implementation of the scheme used for the preparation of arbitrary states of polar- 
ization ququarts is illustrated in Fig. [TJ This scheme employs two nonlinear crystals cut for non-collinear 




Figure 1: A scheme for preparation of an arbitrary pure ququart state using two type-I crystals. 

frequency non-degenerate type-I phase-matching when both (signal and idler) photons are equally polar- 
ized. CW laser is used as a pump, the direction of linear polarization of the pump radiation is specified 
using a half- wave plate (WP), and the relative phase between horizontal and vertical components is 
controlled by tilting a pair of quartz plates (QPs). The state of biphotons generated in the process of 
spontaneous parametric down-conversion under these conditions has the form 

\^) = ^JI\H 1 )\H 2 ) + ^T~il\V 1 )\V 2 ) 1 (3) 

where indices 1 and 2 correspond to different frequency modes. The state of each individual photon in 
the pair is determined by the density matrix 



p J =^\H j ){H J \ + [l-^Vj)^, 



(4) 
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where index j — 1 and 2 numbers the subsystems. 

In this way, we can prepare a state with the required values of the coefficients in the Schmidt de- 
composition but in the fixed (H-V) basis. Transformation to an arbitrary basis |(2J), however, presents 
no special problems. Actually, an arbitrary polarization state of the qubit, i.e., the polarization state 
of a photon in a single field mode, can be obtained from a specified state by means of transformations 
performed with the use of consequently placed quarter- wave plate (QWP) and half- wave plate (HWP). 

Let this transformation be Uj. Then, we have \Hj) — ^> \Aj) and simultaneously \V ) -A \Bj). The 
phase 4> can be eliminated by means of plates QP. As a result, we obtain the transformation 

W) \Ai) \A 2 ) + ^—Ji IB,) \B 2 ) , (5) 

which ends the procedure of preparing an arbitrary state of a ququart. For experimental convenience, 
modes 1 and 2 are spatially separated due to the non-collinear phase-matching. After the transformations 
in eac h sp atial-frequency mode are performed, two modes are combined on a dichroic beam splitter 
(DBS)[31| to obtain a single spatial mode. 

The use of frequency non-degenerate regime is an essential feature of this scheme, making it possible 
to form a collinear beam of biphotons without losses. This is really convenient for practical applications 
in various quantum information protocols. However, in practice, the preparation of entangled states of 
ququarts in the scheme with two consequent crystals involves certain problems. These problems are 
associated with frequency dispersion effects leading to distinguishability of the photon pairs generated in 
the first and second crystals and, correspondingly, to the decrease in purity of the prepared state. 



III. EFFECTS OF FREQUENCY DISPERSION AND THEIR COMPENSATION 



Only two frequency modes are singled out in expression ((3]). Experimentally, this corresponds to the 
filtration of the initial spectrum with ultimately narrow filters. When using filters with a finite bandwidth, 
it is necessary to take into account the structure of SPDC frequency spectrum. This leads to a more 
complex expression for the state vector of the biphoton pair. Let the axes of the first and second crystals 
be oriented vertically and horizontally, respectively. Let us consider the case of continuous pumping that 
is linearly polarized at 45°. In the plane wave pump approximation, the individual states of biphoton 
pairs generated in each crystal of length L can be written in the form (23| 



|Vi,a) = J dndOxdOze 1 ^ sine I o^^.wi + n)4- |H (0 a ,«2 - fi) |vac) , (6) 

where A z (f2) = fc (cJi + 0,6*1) + k a (ui2 — 0,6*i) + k e {uj p ) is the longitudinal phase-matching detuning 
(the pump is an extraordinary wave, and the signal and idler photons are ordinary waves in the crystal) . 
We will analyze the case in which the angular modes are fixed; i.e., they are selected with the use of 
narrow apertures. Under these conditions, integration over angles is not performed; hence, we can restrict 
our analysis to the examination of the frequency spectrum. If the pump coherence length exceeds the 
distance between the crystals [Uj], the corresponding amplitudes add coherently and the total state is the 
superposition 



- y 

V %={ff,V} 



J dQF{n, p)al(oJx + 0)o+(w 2 - fi) |vac) , (7) 



where p numbers the polarization states, 

F(Q, H) ee F(Cl) = sine 1^^) (8) 

and 

F(Q,V) =F(n)e iv( -^ (9) 

are the biphoton amplitudes in the first and second crystals, respectively. 

The presence of additional relative phase ip(fl) is determined by two factors: the phase difference 
between the polarization components of the pump appearing after the first crystal and the appearance of 
an additional phase for the photons generated in the first crystal during their passage through the second 
crystal 

<p(Sl) = {k e (u p ) - k Q (uj p ))L + (fc e (wi + Q) + k e {uj 2 - Q))L. (10) 
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The first term on the righthand side of expression (|10[) does not depend on the phase-matching frequency 
detuning f2 and can be easily compensated by introducing the corresponding phase shift between the 
polarization components of the pump. The second term is related to the frequency dispersion and 
determines the appearance of different phase shifts for different frequencies within the spectral linewidth 
of SPDC. Physically, this means that if a fairly wide part of SPDC spectrum is selected, its individual 
components will be added with different phases, which is equivalent to loss of coherence and to generation 
of a mixed ququart state. The polarization-frequency state will exist only in a Hilbert space of a higher 
dimension D > 4. This dimension can be estimated, for example, by using the Fedorov parameter [24j 
adapted to polarization-frequency distributions; however, this problem is beyond the scope of this work. 

Let us consider the influence of dispersion effects on the properties of generated polarization state. The 
polarization density matrix, which corresponds to the state |^>), can be obtained by taking the partial 
trace over the frequency variables, 



Ppoi = Tro \tp) (tp\ 



- E 

2 ^ 

p={H,V} 



cMF(n, P )F*(n, P ) \pp)( PP \ 



(11) 



in the obvious notation \pp) = a£(wi + ^)ap( w 2 — ^) |vac). We obtain the following expression for the 
polarization density matrix: 



Ppol = <FF| + \VV) (VV\ 



dn \F(n)\ 2 e -'V( n ) |ff) (vv\ + J dn \F(n)\ 2 \vv) (ff|), 



(12) 



which, generally speaking, corresponds to a mixed polarization state. The pure state can be obtained 
only in the case where tp(fl) = const. 

To the first order in frequency detuning 51, we have 



A,(n) 



dk (u) 



( dk (uj] 



n = 



1 



(w 2 ) 



o = c n, 



(13) 



v>(n) = 



dk e (io) 
dui 



dk e (u] 
dui 



m = 



'(wi) vf{w 2 ) 



m = c e m. 



(14) 



Here, v^ r e (uj) = ^ gf" ^ is the group velocity of wave packets with the corresponding central frequencies. 

Expression (| 13[) determines the shape of SPDC spectral line: |F(f2)| 2 = sine 2 (C Lil/2). In our case, 
it is more convenient to consider the temporal representation and to turn from spectral characteristics 
to temporal correlations. Second-order correlation function and spectrum are related by the following 
expression (26| : 



G (2) (r) = 



dQF(n) cos(flr) 



(15) 



When the function F(il) is symmetric, we obtain a conventional Fourier transform. Therefore, in our 
case, the correlation function of biphotons from each crystal has the form of a rectangle with a width 
of Ti = C L. The presence of the phase shift (p(£l) in temporal terms implies a relative shift of the 
correlation functions by t 2 ~ C e L (see Fig. [2]). 

The relative shift in the correlation functions of biphotons from the first and second crystals leads to 
their temporal distinguishability, and hence, to a decrease in visibility of polarization interference. Let us 
consider this problem in detail from the experimental point of view. Scheme of the experimental setup 
is depicted in Fig. [3J Helium-cadmium laser with a wavelength of 325 nm was used as a pump. Vertical 
polarization was selected from the initially unpolarized radiation with the use of a Glan prism (V), after 
which it was rotated with a half-wave plate. Relative phase <f> between the polarization components 
of the pump radiation was introduced by two quartz plates with vertically oriented optical axes. The 
effective thickness of these plates was determined by the relative tilt angle 0. The parametric down- 
conversion occurred in two BBO crystals (2 mm thick each) cut for the collinear frequency-degenerate 
phase-matching. The directions corresponding to the operating wavelengths Ai = 600 nm and A2 = 710 
nm were separated by the apertures (A). Then, the beams were combined on a dichroic beam splitter 
(DBS), transmitting radiation at a wavelength of 600 nm and completely reflecting at 710 nm. In the 
measurement part of the scheme, the pairs were separated in frequency with the use of a similar beam 
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splitter. Film polarizers oriented at 45° with respect to the vertical and interference filters with central 
wavelengths of 600 nm and 710 nm and a bandwidth of 10 nm (IF) were inserted in each channel of the 
measurement scheme. The coincidences of photocounts of detectors Dl and D2 were recorded using a 
coincidence circuit (CC) with a window of 2 ns. 

In ideal case, when the polarization state of a biphoton has the form 

|V) = 1/V2(|friff 2 > + e*\ViV 2 )), (16) 
the coincidence counting rate in this scheme is given by the expression 

R c oc | (+45i, +45 2 | (|ifii? 2 ) + \ViV 3 )) f oc 1 + cos 0, (17) 

where the phase <f> changes with variation in the tilt of the quartz plates in the pump beam. In practice, 
since the pairs generated in the first and second crystals are partially distinguishable due to the incomplete 
overlap of the corresponding correlation functions (see Fig. [2]), the visibility of the interference pattern, 
which is defined as 

t r RC max -Rc mm / 10 \ 
~R T~R ' ^ ' 

^C max ' ^C min 

differs from unity substantially. By using (|12[) for the polarization density matrix, we can derive the 
following expression for the coincidence counting rate: 

R(4>) oc 1 + 1 J dQ \F(Q)\ 2 cos(p(n) -<t>) = 

1 + Vcos{4>), 
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where the visibility of the interference pattern (to the first order in detuning) may be written in the form 

r (~< tci 

V = / dOsinc 2 (-^— )cos(C e iO). (20) 

The integral is calculated simply and leads to the following expression for visibility: 

V = l-^ = 1-^-, (21) 

which is consistent with the intuitive expectations illustrated in Fig. O 

For a fixed frequency spectrum, which, in our case, is determined by the bandwidth of the interference 
filters (it appears to be somewhat smaller than the width of the SPDC spectrum, which, in the regime 
used, is approximately 12 nm), the visibility can be increased by introducing an additional birefringent 
compensator QP into the beam with a shorter wavelength. Quartz compensator with vertically oriented 
optical axis introduces a delay for the vertically polarized photon which is represented by the expression 

Tcomp = [( qr , T oT/ \ )Lq\ , (22) 

where v 9 J^ are the group velocities of the ordinary and extraordinary waves in quartz, respectively, and 

L q is the compensator length. By choosing the compensator length to satisfy the condition T comp = T2, it 
is possible to achieve complete overlap of the correlation functions and, consequently, maximum visibility 
of the interference pattern. 

Figure 2] shows the experimental dependence of the visibility of interference pattern on the length of the 
quartz compensator used. The dotted curve bounding the triangle corresponds to the convolution of two 
rectangular correlation functions, i.e., the functions determined by the shape of the SPDC spectral line 
only, without filtration. However, the experiment was performed using interference filters characterized 
by a Gaussian shape of the dependence of transmission coefficient on the wavelength with the width 
comparable to that of the SPDC spectrum. Therefore the correlation functions should have a shape 
somewhat different from rectangular. When the frequency filtration is taken into account, the expression 
for the biphoton amplitude additionally acquires the factor corresponding to the transmission amplitude 
of the interference filters. By choosing this factor to be Gaussian, we obtain 

F(fi) ex exp(- 1 X ^ )exp(- V 2 )sinc > ( 2 3) 

where Aw is the interference filters bandwidth. The solid curve in Fig. [4] represents the numerically 
calculated dependence of the visibility on the group velocity dispersion compensator thickness taking the 
10 nm interference filters into account. 

The calculation of width and shift of the correlation functions in the case under consideration leads to 
the following values: t\ = 96 fs, T2 = 87 fs, which corresponds to a compensator thickness of 2.6 mm, in 
good agreement with experimental data. The small difference between the experimental position of the 
maximum in the visibility (2.3 mm) and the calculated value is most likely explained by the errors in 
the values of group velocities used. The characteristic interference patterns obtained with and without a 
compensator are shown in Fig. [5j 

The the non-absolute visibility even with the use of a compensator of required thickness is explained 
by two factors. The first factor is "technical" and lies in the inaccuracy of selecting the necessary spatial 
modes. The second factor is the influence of the terms quadratic in the frequency detuning (p(Q) in the 
expansion of f(fl), which lead to unequal broadenings of the correlation functions and, consequently, to 
the deviation of visibility from 100% even in the case of ideal compensation for the part that is linear in 
the detuning (25|. Indeed, to the second order in the detuning from the degenerate frequency, we have 



= « Hnr + ttt L ° 2 = B ? m2 - ( 24 ) 




If we denote £ = C a Lfl 11 D = than for the coincidence counting rate in the case of a complete 
compensation of the first-order effects, we have 

R{4>) oc 1 + — (V c cos <j> + V s sin 0) , (25) 

27T 

where V c = J (i£sinc 2 1 cos D£ 2 nV s = J cosine 2 ! sin Dl; 2 . For the crystals used in the experiment, we 
have D = 0.0193 and numerical calculation gives the following value of interference pattern visibility: 

V = V c n 8 = 0.89. (26) 
2n 
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Figure 4: Dependence of polarization interference visibility on thickness of group velocity dispersion compensator. Dashed 
curve is the theoretical dependence not taking filtering into account. Solid curve is numerically calculated dependence 
taking 10 nm bandwidth interference filters in each channel into account. 
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Figure 5: Interference patterns obtained without compensation and with a 2.33 mm thick compensator. 



It should be noted that such a low value [33| of visibility is the lower limit, corresponding to complete 
absence of frequency filtration. The results of numerical calculation of the dependence of visibility on the 
bandwidth of interference filters used in the case of a complete compensation are presented in Fig. [6l It 
can be seen from this figure that, even when sufficiently broadband interference filters are used (with a 
bandwidth of less than 30 nm), the values of visibility can reach 0.95 or higher. 

The low visibility of polarization interference corresponds to poor quality of entangled states prepara- 
tion. The currently obtained results of statistical state reconstruction for the state 

\$+) = l/y/2(\H 1 H 2 ) +\V U V 2 )), (27) 
which corresponds to the maximum of the interference pattern shown in Fig. lead to a fidelity, defined 



20 40 60 80 100 

filter bandwidth (nm) 



Figure 6: Dependence of polarization interference visibility on the interference filter bandwidth in the case of complete 
group delay compensation. 



as 

F = \{^theory\^exp)\ 2 , (28) 

not exceeding 75%. Nonetheless, we do not see, in principle, any restrictions from above on this quantity. 
Higher values of fidelity can be reached by means of careful selection of spatial modes and, possibly, 
with the use of narrower band frequency filters. However, the mere fact that we observe the two-photon 
interference with a nonzero fidelity suggests that the preparation of entangled states according to the 
proposed scheme is realistic. Regarding improvement of the prepared states quality, it is not a crucial 
and rather technical problem, wich we are working at now. 

IV. CONCLUSIONS 

We have considered a method for preparing arbitrary entangled states of polarization ququarts in non- 
collinear frequency non-degenerate type-I SPDC regime according to the scheme with two consequent 
crystals. We have elucidated the role of effects that are related to the frequency dispersion and which 
lead to a decrease in purity of the prepared states and, as a consequence, to a decrease in visibility of the 
two-photon polarization interference. It has been shown that, in the frequency non-degenerate regime, 
the compensation of group velocity dispersion, which was previously used in schemes with pulse pumping 
only, is necessary even in the case of continuous pump with large coherence length. Methods providing 
this compensation have been proposed and demonstrated experimentally. Moreover, in the case of a 
complete compensation, the maximum value of visibility limited by dispersion effects of higher orders in 
the frequency detuning has been estimated. 

We would like to emphasize that the application of frequency non-degenerate biphotons, unlike the 
majority of earlier proposed and currently known schemes (using frequency degenerate biphotons and 
spatial modes selection), offers a number of advantages. The main advantage is the possibility of forming 
states, including entangled ones, which are localized in a single spatial mode. This property is definitely 
desirable from the viewpoint of possible practical applications, when two-particle states must be trans- 
mitted through a single communication channel. Thus, current work has revealed physical restrictions 
on the quality of prepared states and demonstrated methods for its improvement using the particular 
experimental scheme. 

We would like to thank D. Ivanov for checking the calculations. 
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